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e Average CFD FM accuracy was 2.4% (2009) 
¢ It was believed that poor rotor wake resolution was responsible 
« This lead to research in off-body (OB) adaptive mesh refinement (AMR) 
¢ In 2011 (Chaderjian/Buning): CFD FM predicted with 0.2% for V22 TRAM 
* Vortex wake resolution had no effect (10%, 5%, and 2.5% c,,,) 
¢ Rather, it was crucial to 
= Adequately resolve the formation of the blade-tip vortex 
~+ Fine surface mesh near rotor tip and high-order spatial accuracy 
= Maintain a physically realistic turbulent eddy viscosity in the vortex wake 
+ Detached eddy simulation (DES) turbulent length scale 
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Motivation 


In 2012 (Chaderjian/Ahmad): UH-60A rotor in hover and forward flight (C8534) 
¢ Airloads did not depend on rotor wake resolution 
Both studies did not involve significant blade/wake interaction 
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¢ Animportant question remains 


¢ How are the forward-flight CFD airloads affected by rotor-wake 
resolution when there is significant blade/wake interaction? 


“ Practical engineering issue: High resolution wakes are too expensive 
for most engineering applications 


¢ Two examples for a UH-60A rotor in forward flight are examined 
¢« Blade vortex interaction (BVI), flight-test counter C8513 
¢ Dynamic stall with BVI, flight-test counter C9017 

e Also examine 2D dynamic stall 
¢ Discuss similarities and differences in 2D and 3D dynamic stall 


Outline os 


Flight-Test Data 
Numerical Approach 
Numerical Results 
* BVI-—UH-60A (C8513) 
* Dynamic stall 
= 2D example 
= 3D UH-60A (C9017) 
Concluding Remarks 
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Flight-Test Data/CFD Validation 


Joint NASA/U.S. Army UH-60A Airloads Database (1993/1994) 
¢ Airloads at various radial locations along the rotor blade 
Bousman’s qualitative analysis of dynamic stall (AHS Journal/Oct. 1998) 
“« He examined the time history of blade pressures to judge when 
= Moment stall: Formation of dynamic stall vortex at blade leading edge 
= Lift stall: When dynamic stall vortex passes over blade trailing edge 
= Flow separation at blade trailing edge 
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(CFD/CSD Loose Coupling) 


OVERFLOW 2.2L — CFD Flow Solver 
¢ Solves the time-dependent Navier-Stokes equations 
* Structured overset grids 
* 2¢-order dual time accuracy (At=%° rotation, 60 subiterations) 
Loose « At least 2.3 subiteration residual drop for all grids 


eee * 5 order spatial accuracy (central differences/artificial dissipation) 
very 


VarevGlinon * Hybrid RANS/LES turbulence model 

= Spalart-Allmaras one-equation turbulence model 
= DDES length scale 
» SARC rotation/curvature correction 
= Y*<1 at body surfaces 

CAMRAD II — Helicopter Comprehensive Analysis Code 

¢ Provides rotor-blade structural deflections 

¢ Provides trim control angles at the blade root 
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(Near-Body and Off-Body Grids ) 


Cartesian Off-Body Grids Rotor-Blade Grids 


Blade Tip 


Level-1 Cartesian Grid Level-1 AMR 
Cartesian Grids 


Trim Tab 


Level 2 


¢ Rotor blades/Hub use O-mesh topology 
¢ Off-body grids use Cartesian grids with adaptive mesh refinement (AMR) 
* Rotor wake captured only with Level-1 grids (10%, 5%, and 2.5%c,,,.) 
¢ No interpolation throughout the resolved rotor wake of interest 
* Has same resolution and coincident overlaping grid points 
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Flight-Test Data 
Numerical Approach 
Numerical Results 
* BVI —UH-60A (C8513) 
¢« Dynamic stall 
= 2D example 
= 3D UH-60A (C9017) 
Concluding Remarks 


Sy BVI Flight Counter C8534 ke 
NASA’S OVERFLOW Navier-Stokes CFD Code 
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Wake-Grid Spacing 
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Three AMR Wake-Grid Resolutions 
BVI Flight Counter C8534 


NASA’s Pleiades Supercomputer 
5,628 Broadwell CPU Cores 


AS=10% Cy, 


500 Grids 


87 Million Grid Points 
4.6 Hr/Rev 
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297 Million Grid Points 
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12,000 Grids 
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Effect of Wake-Grid Resolution on Airloads ee 
BVI Flight Counter C8534 


¢ Good overall agreement with flight-test data 


¢ OB resolution has very little effect on airloads! 
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Flight-Test Data 
Numerical Approach 
Numerical Results 
* BVI -—UH-60A (C8513) 
¢« Dynamic stall 
= 2D example 
= 3D UH-60A (C9017) 
Concluding Remarks 


CFD Animation of Two-Dimensional Dynamic Stall ee 
= 10° + 10° sin(2kt — 4), k= = 0.1 


How Most People 
Think of Dynamic Stall 
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Vortex passing airfoil trailing edge — (lift stall) 
Reversed flow 

¢« As vortex forms at leading edge 

¢« As vortex passes trailing edge 
Three stall events, each smaller than the previous one 
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T Upstroke 


Two-Dimensional Dynamic Stall C——- 
Pai i eine as 
2V. 
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Force/moment time-history 
indicates 3 stall events 
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2-3 typical 

2D characteristics identified 
experimentally (McCroskey 
et al., 1976) 

Many feel 2D captures the 
essential elements (Tan & 
Carr, 1996) 

It will be shown that 2D does 
miss some important 3D 
dynamic stall characteristics 
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¢« Dynamic stall 
= 2D example 
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Concluding Remarks 
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High-Resolution Dynamic Stall (C9017) 
NASA’S OVERFLOW Navier-Stokes CFD Code 
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There is BVI 
It is affecting the dynamic stall process 
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Three AMR Wake-Grid Resolutions 
BVI Flight Counter C8534 


NASA’s Pleiades Supercomputer 
5,628 Broadwell CPU Cores 


—_. 760 Grids 
83 Million Grid Points 
: 4.5 Hr/Rev 
| eee 3,200 Gri ae - AS=5% Cig 
241 Million Grid Points 
6.2 Hr/Rev 
— ane 


1.3 Billion Grid Points 
28.5 Hr/Rev 
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Effect of Wake-Grid Resolution on Airloads Sa 
Dynamic Stall Flight Counter C8534 
Good overall agreement with flight-test data 


More high-frequency content, but little effect on airloads! 


This suggests AS=10%c,,, adequate for engineering design airloads 
r/R=0.675 r/R=0.865 : r/R=0.965 
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Closeup View of 3D Dynamic Stall With BVI os 
Inboard and outboard separation, with attached flow in between 
3D Vortex rings emitted due to Helmholtz vortex theorem 
¢~ Different from 2D Vortex lift-off 
Vortex path altered due to separation 


* Can effect aeromechanics of following rotor blades 
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Tip vortices from Blades 2 and 3 do 
not disrupt the flow on Blade 1 


Blade 1 
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Nomenclature 
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Closeup View of Dynamic Stall a 


e Flow separates outboard of vortex, 
remains attached inboard of vortex 


¢ Separation moves with the vortices 
e Tip vortices from Blades 2 and 3 
appear to trigger dynamic stall 


Blade 1 at y=225° 


Blade 1 


Separation m 
Blade 2 
Vortex Blade 3 


Vortex 


Rotating Blackhawk Rotor Blade 
But What Happens When a Vortex Passes Over the Blade? 
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fe Lower Incidence 
Higher Incidence 


Velocity Vectors Relative to Rotor Blade ee 


Inboard of Vortex 


Outboard incidence is greater NG 
than inboard incidence by at 
least 10 deg 


¢ This explains why stall occurs 
outboard of the vortex and 
reattaches inboard of the 
vortex 
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Closeup View of a Single Blade er 


(Same Blade Motion and Aeroelastic Deflections) 
¢ No outboard separation in the 3 quadrant! 
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« This confirms vortex-induced dynamic stall 
¢ Inboard separation due to freestream reversed flow 
¢ Separation along entire blade in 4 quadrant, due to blade incidence 
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BVI-Induced Dynamic Stall (C9017) Bee 


2D Wind-Tunnel Experiment 


3.5 chords 


e First observed experimentally for a 2D airfoil 
* 38° European Rotorcraft Forum: Zanotti, 
Gilbertini and Mencarelli 
* Similar explanation of how a vortex 
triggers dynamic stall 


Vortex generator Oscillating airfoil 


First time observed for an actual helicopter rotor 
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al Comparison of CFD With Qualitative Flight-Test Analysis “<-> 
(Dynamic Stall, C9017) 

¢ Polar plot 

* Bousman’s moment stall, lift stall, and trailing-edge 


separation 
—e— EXP Moment Stall 
*« There are two stall events an 
5 x —— CFD Outboard Vortex 
* CFD Outboard andinboard vortices 9 t&= araikeceiveces |" " 
1 
= Outboard vortex initially moves inboard then ae: 
outboard 


* Tracks stall closely up to 270°, where it drops below 
the blade and has little influence 


* Inboard vortex only moves inboard 
e Flight test does not indicate inboard reversed flow 


¢ It must be there, but loads are light and pressure 
data sparse (Bousman) 


Vorticity Magnitude 
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Blade Reference Frame 


Blade Stalls As It Retreats 
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Conclusions 


Good overall comparison between CFD airloads and flight-test measurements 
for BVI and dynamic stall cases 


¢ wake grid resolutions were AS=10%, 5%, and 2.5% Cin 


Refining rotor wakes beyond engineering resolution (AS=10% C,,,) did not 
significantly affect the predicted airloads, even with blade/wake interaction 
« This suggests that airloads engineers may use the coarser wake-grid 
resolution (AS=10%c,;,) for hover and forward fight simulations provided 
= The CFD tip-vortex is accurately formed using a combination of fine 
surface mesh at the blade tip and high-order spatial accuracy 
« Use of a hybrid RANS/DDES turbulence model 
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Conclusions ——— 
(Continued) 


¢« Differences between 2D and 3D dynamic stall 
= 3D vortex rings are emitted rather than a simple 2D leading-edge vortex 
= Dynamic stall flow separation can alter the path of a BVI vortex 
= Vortices passing over the rotor blade caused BVI which triggered dynamic stall 
+ This phenomenon has been observed in a 2D wind-tunnel experiment 
= Mechanism for BVI-triggered dynamic stall 


+ Induced velocity field by other blade-tip vortices changed the relative angle of 
attack of the stalling rotor blade 


o The blade AOA increased outboard of the BVI vortex, causing flow separation 
o The blade AOA decreased inboard of the BVI vortex, resulting in attached flow 
* The successful modeling of 3D dynamic stall with BVI should include an accurate 
prediction of blade-tip vortex trajectories 
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